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1 Introduction 



A paper by E. T. Tomboulis that was posted last year on the arXiv pQ claims 
to present a rigorous proof of quark confinement [2\ in 4D lattice gauge the- 
ory. This problem has been fascinating many physicists since the last century 
and still remains an important open problem; in fact, its resolution probably 
would represent a first step towards the solution of one of the millennium 
problems posed by the Clay Mathematics Institute [3] . In a comment posted 
on the arXiv [4J we pointed out that in our opinion the proof in [1] depended 
on a claim which remains to be proved. The purpose of this public comment 
was to prevent confusion in the community about the question whether con- 
finement on the lattice was proven or not. Tomboulis then posted a reply [5] 
on the arXiv stating that the presumed gap found by us was simply due to 
a misunderstanding of his argument. 

We appreciate that Tomboulis gave a concise and clear exposition of the logic 
of his argument in his reply. But we still think that our original criticism 
stands and in this note we try to give a more detailed explanation of the main 
open problems we see in Tomboulis' approach. We also discuss possible ways 
to complete his proof. 

2 Tomboulis' claimed theorem 

To make this note self-contained, we summarize Tomboulis' notation and his 
arguments with some simplifications: 

1. A C Z 4 is a box in Z 4 of size L\ X . . . x L4 with center at the origin, 
where Lj = ab ni , >> 1, i = 1, ... ,4 and a, b are positive integers 
larger than 1. 

2. C Z 4 is the box in Z 4 of each side length Lib~ k obtained from A 
by k steps of the renormalization transformation of Migdal-Kadanoff 
type (A = A^). Periodic boundary conditions are employed in all 
directions of /1 = 1, ... ,4. (actually periodic boundary conditions in 
the directions of fi = 3,4 would be enough.) 

3. f({cj},U) = 1 + Ylj^o c jdjXj(U) where dj is the dimension of the 
representation Xji an d we assume U G G = SU(2). Moreover U = 



2 



Up — Ylbed-p Ub for a plaquette p C A which consists of bonds b G dp 
oriented positively with respect to p. 



We start with 



exp 



fxi/ 2 (£/) 



F o (0)f({ Cj },U) 



f({ Cj },u) = i + Y,cAxj(u) 



^0(0) 



exp 



;Xi/2(U) 



dU 



(2.1a) 
(2.1b) 

(2.1c) 



where dU is the Haar measure of G — SU(2). We then apply the renor- 
malization group (RG) recursion formulas of Migdal-Kadanoff (MK) type 
[6] with some modifying parameters. If there are no such parameters (the 
standard recursion formula), we have 

f^\U) = f({ Cj (n-l)},U) 
Mo 



f (n) (U) = f({ Cj (n)},U) 



where 
f {n \U) 

F (n) 
or equivalently 



FJn 



f^iuu^f^iu^u,) ■ • ./ (B_1) (tfa) fY[dU k 

S v ' 

b 2 

(2.2a) 

[f(n-i)(u}f dU \ (2.2b) 



Fn(n) 



Fj(n) 



[f (n-D {u) f 



Xj(U) 



d; 



dU 



b 2 



(2.3a) 
(2.3b) 



in terms of the coefficients of the character expansions. Then [7] : 
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Theorem 2.1 For D < 4 and for G = SU(N) or G = U(N), 

lim Cj(n) = for j ^ 

n^oo 

These recursions are just approximate and yield upper bounds for the parti- 
tion functions [Ij. Tomboulis [I] then introduces two parameters a G (0,1] 
and t > and a function /i(a, £) G (0, 1], designed so that this transformation 
becomes numerically exact for the partition function: 



Z = ! dU A Y[f({ Cj },U p ) (2.4a) 

J pCA 

= [F (1)^]I A(1) I fdU AW J] fm(a)},U p ) (2.4b) 
J pcAC 1 ) 

= Zt(c{a),t) (2.4c) 
where dC/ A = H beA dU b , 

Z x {c{a),t) = [F (l) h ^ AWl Zi (2.5a) 
Z x = [dU AW H /({£,•(«)}, E7 P ) (2.5b) 

and as usual, {7 P = ribedp ^ b are plaquette actions defined as the product of 
group elements £/& = Ur X)X+e \ G G attached to the (oriented) bonds 6 G A. 
Moreover 

h(a,t) = exp[— £(1 — en) /a] (2.6a) 
Cj{a) = cf ] (a) = acj{l) (2.6b) 

«,<!) = gg, (2.60 

(where without loss of generality we have chosen one particular form of the 
function h, as suggested in []]). Given £ > 0, a G [0, 1] is chosen so that the 
above relation becomes exact; thus a has to be considered as a function of £. 

The author of pQ then introduces a vortex sheet V — {v C A} [SI [9] which is 
a collection of plaquettes {v = {x + n 3 e 3 + n 4 e 4 , x + e x + n 3 e 3 + n 4 e 4 , x + 
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e 2 + n 3 e 3 + n 4 e 4 , x + e 1 + e 2 + n 3 e 3 + n 4 e 4 }, rij = 0, 1, , . . . , L u (i = 3, 4)}, i.e. 
a plaquette p = (x , x + ei, x + e\ + e 2 , x + e 2 ) in an a?i — x 2 plane and its 
translates along the axis normal to p (say, 3rd and 4th axis). Following p] 
we define 

= JdU A Y[f({ Cj },(-iyMu p ) (2.7) 



pCA 



where 



and then 



»0) = \\ZVv (2 ' 8 » 



= fdU A ]J (l + J2cjd jXj (U P )) 

J P CA\V j^O 

xII^+E^ViXi^)) (2-9) 

qCV j^O 

Note that the position of the vortex sheet V C A can be freely moved in the 
X\ — x 2 plane by gauge invariance. 

First of all, note that the coefficients of XjiP) in Yl (1 + J2j^o c jdjXj(.U p )) 
and 

n(i + 5>dix*(Op)) + ri( i + e^^- 1 )^^)) 

p jyo p jyo 

are nonnegative. The main claims in [1] are: 

Claim 2.1 ('ij There exist t, t + > and functions a and a + such that 

z { a\{cA) _ 47n)({gS w) (Q + (* + ))» 

^({c,}) Z AW ({4 n) (a(t))}) 



c 



j 



a(t)) = a{t)cj{n). 



(2) For large n such that {cj(n)} are sufficiently small, there exists a > 
(the same for numerator and denominator!) such that 



ZaUc,}) " Z A(n) ({cf\a(Q)}) 
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If this were correct, the following would be true: 




tends to the high temperature fixed point (i.e. {cj(n)} — > 0) as n — > oo if the 
dimension is < 4, whether G is abelian or non-abelian (see the remark below 
and [7]), this would mean strict positivity of 't Hooft's string tension and thus 
establish permanent confinement of quarks in the sense of Wilson at least for 
all values of the bare coupling constant [9j [10] in 4 dimensional lattice gauge 
theory, thereby solving part of a longstanding problem in modern physics 
(the missing part, as Tomboulis makes clear, would still be the proof that 
there is a continuum limit in which confinement persists). 

As pointed out in our earlier note [3], this cannot possibly be correct, however, 
since there exists a deconfining transition in 4D lattice gauge theory based 
on abelian gauge groups. But it is meaningful for the future study of the 
confinement problem to ask what is going wrong in reaching this conclusion. 

The original idea to study vortex free energies in order to understand con- 
finement was formulated by Mack and Petkova in [9]. They fix boundary 
variables Gqa = {gt G G;b — (x, x + e M ) G <9A} and take the maximum of 
Z A ^/Za over Gqa- On the other hand 't Hooft's string tension is defined 
using (twisted) periodic boundary conditions at <9A, see [10J. The following 
are well-known facts: 

Theorem 2.2 Define 



(1) LetG = U{\) orG = SU{2). Then there is a fa > such that for (3 < fa 



(area decay law), a > is, if chosen maximally, is 't Hooft's string tension. 
(2) Let G — U(l). Then there is a (3q < oo such that for (3 > (3q 



F = - log 




Z A 




F < const exp[— CTL1L2] 




F > const 



1 



(deconfinement) . 
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(1) is a standard result of the convergent high-temperature expansion (see for 
instance [33]). The analogous statement for SU(N) holds with Z^ replaced 
by Z w , the partition function twisted by an element to of the center of SU(N). 

(2) is a variation of the results of [11] and [12] (see [T3]). 

Remarks 2.1 (1) The introduction of < a < 1 into 1 + ^2j^ Q CjdjXj{U) 
does not violate conditions (positivity, analyticity, class functions etc.) on 
f( n \v) in l2H since 

1 + aCjd jXj (U) = (1 - a) + a f 1 + J2 c J d jXj(U) J 

and < a < 1. JTien {cj(n) > 0} tends to as n — > oo (if r = \, see Section 
3.3). 

(2) The conjecture raised in JM/ is proved rigorously in fiTj]/ , Namely 'tHooft's 
string tension is smaller than or equal to Wilson's string tension. 

(3) The center of G = U(l) is again U(l). Then we consider U(p) = cos(# p ) 
and Uu{p) = cos(9 p + uj) for p G V . U^ip) — —U (p) for lo — it. 

(4) The argument in [11] uses the Fourier transformation of e 13 ^ 9 and 
Poincare's lemma. Since A is a torus, Poincare's lemma has to be replaced 
by the Hodge decomposition in order to adapt the arguments in f77]/ or IIB] . 
See Hf and \T^. 



3 Tomboulis' Proof revisited 

We now analyze the arguments in [lj in more detail. First, using the fact 
that the partition function Z = Z\ increases by the Migdal-Kadanoff (MK) 
recursion formula [T], Tomboulis introduces two interpolation parameters a 
[Z increases as a /* 1) and t (the factor [F (n)] /l, ^ Q '^ decreases as t increases). 
Then he claims that there exist functions a(t) and a + (t) such that 

Za({cj}) = [FoWl^^l^l^^CaCt))}) (3.1a) 

Zl({c 3 }) = [FoWl^^l^lZ+^^Ca+W)}) (3.1b) 

where 

z+ = l -{z + z^) 
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and the right hand sides are manifestly independent of t. 

The author of fl] then claims is that there exists a t* > such that 

a(U) = a + (Q (3.2) 

which yields 

zj{{c 3 }) = zum) 

Z A {{cj}) Z AW ({dj}) 

where 

The procedure is then iterated to produce an analogous relation involving 
the n— th iteration of the MK RG transformation. 

Two arguments are given in [Tj for this proof. 

large (3 region: p.24-p.26, based on the arguments in pages 14,15,22, 23 and 
Appendix B. 

small (3 region: p.26-p.27 and Appendix C. 

The first argument is used to drive the system to the high-temperature (small 
(3) region by the MK recursion formulas, and the second argument depends 
on the implicit function theorem to derive the main claim which yields strict 
positivity of 'tHooft's string tension. 

3.1 Outline of the Proof in Jl] 

To understand the proof in in PQ, we streamline the arguments in [TJ by 
extracting key parts from [1] (some equations are added below by the present 
authors). But before doing so, we have to stress that Tomboulis is modifying 
the MK recursion by introducing a further parameter r e (0, 1] (Eq. (2.19) in 
PQ), which is the source of some serious problems, because any choice r < 1 
it changes the MK RG flow for large (3 drastically (see below). 

r is introduced by replacing (I2.3al) with Cj(n) = Cj(n) b2r (this amounts to 
replacing the b 2 convolutions in Eq. fl2.2al) by b 2 r convolutions). 
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Quotation (I): [page 14] 



dt 

where 



«K7(i,»-),i,r) (3.25) 



v(a,t,r) = - a ., (3.26) 

dh(a,t) 

— — h A A(m )(a,rJ 

with <9/i(a:, t)/<9t = —[(1 — a)/a]h(a, t) < and 
with 

i a 

A AM (a,r) = p^^-lo g Z AW ({c J ( ffll « 1 r)})>0 

(3.27) 

[ page 15] We require that 

5' <a Ath (t,r) < 1-5 (3.30) 

. . . by choosing r to vary, if necessary, away from unity in the 
domain 

l>r>l-e (3.31) 

where < e « 1 with e independent of |A (m) |. With (3.30) in 
place, (3.25) and (3.29) imply (Appendix B) that 

^a Ath (t,r)>rn(6)>0 (3.32) 

[page 23] by letting the parameter r vary, if necessary, in (3.31). 
From this, it follows that the analog of (3.32): 

^ t al h (t,r)> V t(5 + )>0,... (4.25) 
holds . . . 
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(italics introduced by us for emphasis). 

Quotation (II): [page 25] ... (5.4), taken at general r, implies 
that . . . 

H, h (t,r)-a A , h (t,r)\<o(J^J (5.6) 

. . .. For any given tf, choose t\ in Z A ^ h {j3, (XA,h(ti), *i) so that 
h(a Ath (t 1 ),t 1 )=h(aX >h (tt),tt) (5.7) 
This is clearly always possible by (3.32) and (4.25), and by (5.6). 

[top of page 26] We may now iterate this procedure (n — 1) deci- 
mation steps . . ., at each step choosing t m and i+ such that 

M«£2(*m),tm)=M«AS B) (0,0 (5-9) 

Carrying out . . . one obtains 

Za J Z A(n) ((3,h,a^(t),t) 

Quotation (III): [middle of page 26] Next, consider (5.10) rewrit- 
ten as 

z * ) Uj-d(^^«S(*),*)/ v ••)}), 

(5.13) 

.... It is then natural to ask whether there exists a value t = 
such that 

Z+ {n _ 1) (P,h,a^ h (t),t) = Z+ (n _ 1) (5.14) 

Note that the graphs of ot^\(t) and a^\t) must intersect at t^\. 
A unique solution to (5.14) indeed exists as shown in Appendix 
C provided 

A A{n) (a,r)>A+ (n) (a,r) (5.15) 
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(The above assertion is to be interpreted as affirmative.) 

[near top of page 27] V.l Consider n successive decimation steps 
performed according the scheme (4.28). Assume that there is 
an n such that the upper bound coefficients cV(n) become suffi- 
ciently small for n > uq. 
[Appendix C, page 44] With the notation 

= i n i^H|AW| lnZ ^ ({5j(n,a)}) (c ' 5) 

we now define 

tt(A,*) = h(a(t),t) + (1 - A)$+„(a+(tj)) + \$+ n (a h (t)) - 

(C.6) 

and consider the equation 

*(A,t) = (C.7) 

At A = 0, eq.(C7) is solved by setting t = t . . .. By the implicit 
function theorem, if grad \I> is continuous and d^f/dt ^ 0, there 
exists a branch t(X) through (0, to) on a . . .. Thus, from (C.8), 
£(A) . . .extends to the solution £(1) > t . 

Quotation (III) is the main result and quotation (II) is used as the bridge to 
reach the main conclusion. The summary of the argument in [1] is: 

(I) If r < 1 - e then a < 1 - 5 < 1 and a + < 1 - 5 < 1. (5 > 0). Thus 

v = da(t,r)/dt and v + = da + (t,r)/dt are bounded from below by 
strictly positive constants r/i > and r] 2 > 0. 

(II) Since \a + (t + ) — a(t)\ = 0(|Ai| -1 ), the existence of t\ such that 

/&(o!A,h(ti),ti) = h i. a l,h^i)^i) 

follows by a shift of t. Moreover it says that one can find t = t m suc- 
cessively to obtain Z+ m) (a+( m ) {t + )) / Z k(m) (a^ (t)). This proves Claim 

O(i). 

(III) After some steps, {cj(n) > 0} are sufficiently small for j ^ and 
we can find t = t n such that aS n \t) = a +( - n \t) by using the implicit 
function theorem. This leads us to Claim [27T1 (2). 
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3.2 Problematic parts in the proof 

The previous proof may look perfect, but there must be something missing 
since this proof seems to work even for G = U(l). We mention some prob- 
lematic points: 

Problem 1. The existence of t m satisfying (5.9) in [1] or in (II) is not 
established. 

Problem 2. The existence of i* satisfying 

= (3-4) 

is not established. 

Problem 3. The implicit function theorem is used to obtain t(X) which 
satisfies ^(A, t) = 0. It is not proven that t(X) can be extended from 
t = to t = 1 even if V t ^ 0. 

First we remark that all these claims depend on the hypothesis that both 
do>( n \t)/dt and da +<yn \t)/dt are bounded from below by strictly positive 
constants 771 > and i] 2 > respectively, uniformly in t G R. But this 
cannot be correct without assuming ct(t) < 1 — 5 (5 > 0) for all t uniformly 
in n or A*™) because of the form of h(a, t) = exp[— t(l — a) /a]. See (I). It is, 
however, very plausible that a is very close to 1 since the original (a = 1) 
MK approximate formulas are very accurate. To ensure those lower bounds, 
the parameter r is introduced. We suspect that this is the origin of the 
contradictions as we explain in the next subsection. 

Secondly we need to find t* such that a(i*) = a + (t*), otherwise we compare 
two different physical systems, which is absolutely meaningless. To do so, 
the author of [lj appeals to the implicit function theorem. But his claim does 
not follow from this theorem directly. We need that the map is contractive 
to ensure that this claim is correct. 

3.3 The origin of the problems. 

The author of [lj tries to ensure that da/dt = v and da + /dt = v + are 
strictly bounded away from zero, and this is the reason for the introduction 
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of the parameter r. The derivatives v and v + contain 1 — a in the numerator 
and A A ( m ) in the denominator which may make v and v + small. Here A A ( m ) 
is (almost) equivalent to the derivative of the free energy by the (inverse) 
temperature and is similar to the specific heat, see (3.27) in Quotation (I) 
or in p]. It is easily seen that A A ( m ) > is bounded from above since 
f({cj}, U) > 0. So the denominator is not dangerous. 

Thus we have to worry about 1 — a — > 0. Since the MK recursion formulas 
are very accurate for large (3, a seems to be very close to 1 for large (3. (This 
is so even though the MK recursion formulas fail to describe the deconfining 
transition of the U(l) model even for a = r = 1.) For this reason, the 
parameter r e (0, 1] is introduced by eq.(2.19) in [lj so that a < 1 — 5. But 
the price to pay for this is that Cj(n) may not converge to (the strong 
coupling fixed point) as n —>■ oo ! 

As remarked above, Cj(n) = dj(n) b2r means that we replace b 2 convolutions in 
(j2.2aj) by b 2 r convolutions. If r is chosen less than 1, the effective dimension 
of the MK recursion formula increases from D = 4 to D > 4: 

2 loer 

D = 4 ^ D = A 



2 log b + log r 



see 016]. The choice r = 1 corresponds to the original MK formula and D = 
4 becomes the critical dimension in the sense that the number of convolutions 
(i.e. b 2 ) is equal to the number of block plaquettes to be glued together (i.e. 
b 2 ). 

If r < 1, the recursion formula drives the system to the weak coupling fixed 
point if (3 is large and a is close to 1. Namely if r < 1 and f3 is large, Theorem 
12.11 does not hold and lim n ^oo Cj in) ^ in general. 

The following pedagogical example makes this clear: Set fo(U) = exp[— j39 2 /2], 
where (3 > and 9 G R. This corresponds to (non-compact) abelian lattice 
QED. Then for the choice of r 6 (0, 1], from f)2.2ap . we get 



fo -> fn(9) = exp 



2r n 



which for r < 1 converges to the weak coupling fixed point (5(9) with suit- 
able normalization). If r < 1, the same result will hold for compact U(N) 
and SU(N) lattice gauge theories for large (3. This is plausible because in 
this regime the Gaussian approximation is good; it has been confirmed also 
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numerically by us and [15J. As an illustration we show in the figure the re- 



sults of numerical iteration for G = SU(2) r = 0.9 and a = 1. There is a 
'critical point' at (3 ~ 4.8: for f3 > 4.8 the modified MK flow goes to the 
weak coupling fixed point, whereas for j3 < 4.79 it still flows to the strong 
coupling fixed point. 




Figure 1: Evolution of Cj/co under Tomboulis' modified MK RG with r = 0.9. 
(3 = 4.80 (left plot), (3 = 4.79 (right plot); lines drawn to guide the eye. 

It should not be difficult to prove this rigorously. These examples show that 
one has to be very careful if r is chosen < 1: depending on (3, r has to be 
chosen close to 1, which in turn will mean that a is close to 1. 

The author of pQ says that "choose r < 1 if necessary ..." in several places 
(e.g. 3rd, 16th lines on page 15, 8th line on page 18, 15th line on page 23, 
etc.) Though it is said that "r < 1 is actually irrelevant for G = SU(2)" in 
the 9th line on page 19, r is chosen < 1 throughout the paper. On the other 
hand, to ensure lim^oo Cj(n) =0, 1 — r will have to be chosen very small 
for large (3. But then the required upper bound a(r,t) < 1 — 5 becomes very 
subtle. 

So we have two choices: 

1. r < 1. The larger (3, the closer to one r will have to be chosen to make 
sure that the flow goes to the strong coupling fixed point. But then it 
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will be very difficult to ensure that a is sufficiently far from 1. 

2. r = 1. All parameters depend on (3 and |A| in a very subtle way in this 
case. (But this is sure to fail for G = U(N).) 

These subtle points are not addressed in pQ. 

Finally there is a problem with proving the existence of a t such that a(t) = 
a + (t) after sufficiently many iterations of the MK recursion formulas or when 
{cj} are small. To obtain the function t(X) satisfying \I/(A, t) = with 
t(0) = 0, the author of [1J appeals to the implicit function theorem. Since 



by the parametrization invariance ( t-invariance) of the partition function 
(this is the definition of a + ), we have 



We assume that the equation \I/(A = 0, t) — is solved by t = t , and the 
equation \I/(A,£) = is our required equation, and we want to know if the 
solution t = t(X) with t(0) = to can be continued to t(l). We have 



logZ+({ Cj }) 



log [FoCij^^i^VC^W)}) 

log F (l)^ + ^)'*)l A(1) IZ 1 + ({5(a + (t / ))}) 



$(A = 0,t) = h(a(t),t) - /i(a + (t/),t 7 ) 



(3.5a) 




({c>(t))})] -logZ + ({c,})) (3.5b) 




(3.6) 



F(s,t(s)) 



(3.7) 



where 




(3.8a) 



1 



|A (1) |logF (l) 

x (logZ+({9(a(t))}) -logZ + ({g,(a 



+ 



(*)})) (3.8b) 
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and 

I — a t 

h t (a,t) — h(a,t), h a (a,t) — —h(a,t) (3.9) 

a cr 

The integral equation t(X) = </>(£) (A) can be solved analytically by the iter- 
ations if F(s,t) is bounded in the region. As Tomboulis pointed out in pQ, 
A > A + if {cj > 0} are small and the high-temperature expansion converges. 
But the condition ^(A, t) 7^ does not guarantee that t(X) can be defined 
for all A. The following example is given by Kanazawa [IB] : 

#(A,t) = e~ t -l + 2X 

To justify the claim in [1], we need contractivity of the map <fi (in a suitable 
set of continuous functions) and to prove it is again not trivial at all. 

4 Discussion 

To sum up, the paper [1] contains several problematic points which remain to 
be proved or confirmed. It seems that the remaining problems are, however, 
not easy to solve. 

If the conventional wisdom of quark confinement in 4D non-abelian lattice 
gauge theory is correct, the alleged theorem in [1] is certainly very plausible 
and may hold for G = SU (N). But it is again a very subtle problem to show 
the existence of t* such that a(t») = a + (t*) for a value of r close enough to 
1 to ensure convergence of the MK RG to the strong coupling fixed point. 

The case of G = U(l) highlights this subtlety: For large (3 and r = 1 such 
a t* does not exist, whereas for r < 1 may exist, but does not imply 
confinement. We think in this note we made clear why it is so subtle to 
establish the existence of t*. 

Though the MK RG recursion formulae cannot distinguish non-abelian groups 
from abelian ones, the velocities of the convergences of {cj(n)}°^ 1//2 to as 
n — ► oo are very different. We are skeptical about the idea that the prob- 
lem of quark confinement can be solved by soft analysis like this, but if the, 
MK RG formulas should play a role in a rigorous proof of quark confine- 
ment in lattice gauge theory, these different velocities would certainly have 
to come into play. Presumably the dependence of r and hence a on /3 and A 
must be clarified, but unfortunately it is very difficult to find out what the 
relationship is. 
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